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ABSTRACT 

We present results from mid-IR spectroscopic observations of two young super- 
nova remnants (SNRs) in the Large Magellanic Cloud (LMC) done with the Spitzer 
Space Telescope. We imaged SNRs B0509-67.5 and B0519-69.0 with Spitzer in 2005, 
and follow-up spectroscopy presented here confirms the presence of warm, shock 
heated dust, with no lines present in the spectrum. We use model fits to Spitzer IRS 
data to estimate the density of the postshock gas. Both remnants show asymmetries in 
the infrared images, and we interpret bright spots as places where the forward shock 
is running into material that is several times denser than elsewhere. The densities we 
infer for these objects depend on the grain composition assumed, and we explore the 
effects of differing grain porosity on the model fits. We also analyze archival XMM- 
Newton RGS spectroscopic data, where both SNRs show strong lines of both Fe and 
Si, coming from ejecta, as well as strong O lines, which may come from ejecta or 
shocked ambient medium. We use model fits to IRS spectra to predict X-ray O line 
strengths for various grain models and values of the shock compression ratio. For 
0509-67.5, we find that compact (solid) grain models require nearly all O lines in X- 
ray spectra to originate in reverse- shocked ejecta. Porous dust grains would lower the 
strength of ejecta lines relative to those arising in the shocked ambient medium. In 
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0519-69.0, we find significant evidence for a higher than standard compression ratio 
of 12, implying efficient cosmic-ray acceleration by the blast wave. A compact grain 
model is favored over porous grain models. We find that the dust-to-gas mass ratio of 
the ambient medium is significantly lower than what is expected in the ISM. 

Subject headings: interstellar medium: dust — supernova remnants — Magellanic 
Clouds 



1. Introduction 



Supernova remnants (SNRs) provide a laboratory to study various aspects of interstellar 
medium (ISM) evolution across the whole electromagnetic spectrum. With the advent of high 
spatial resolution telescopes in both the X-ray and infrared (IR) regimes, it is possible to probe the 
interaction of the rapidly moving shock wave with the dust and gas of the surrounding ambient 
medium and to study the ejecta products of the SN itself. The expanding Shockwave sweeps up 
and heats gas to 10^ — 10^ K, causing it to shine brightly in X-rays. Dust grains embedded in 
the hot, shocked plasma are coUisionally heated, causing them to radiate at IR wavelengths, and 
slowly destroying them in the process via sputtering. Because the physical processes behind X-ray 
and IR emission are related, a combined approach to studying SNRs using both energy ranges can 
reveal more information than either could on its own. 

To better characterize IR emission from SNRs, we conducted an imaging survey with the 
Spitzer Space Telescope of ~ 40 known SNRs in the Large and Small Magellanic Clouds. The 
Clouds were chosen because of their known distance and relatively low Galactic IR background. 
Subsequently, we obtained IR spectra of several of the SNRs in the LMC using the Infrared Spec- 
trograph (IRS) on Spitzer. We report here on IRS observations of two of these, SNRs B0509- 
67.5 (hereafter 0509) and BOS 19-69.0 (hereafter 0519). Both are remnants of thermonuclear SNe 
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Hughes et al.ll 19951) and hav e fast, non-radiative shocks (several thousand km 



Ghavamian et al.ll2007|) . There is no evidence in either for slower, radia- 
tive shocks. They are both located in the LM C, at the known d istance of ~ 50 kpc. In addition, 
both have ages determined from light echoes (IRest et al.ll2005b . with 0509 being 400 ±120 and 
0519 being 600 ±200 years old. They are nearly identical in size, having an angular diameter of 
~ 30", which corresponds to a physical diameter of ~ 7.3 pc. 



In 



Borkowski et al.l (|2006|) hereafter Paper I, we used 24 and 70 /im imaging detections of 



both SNRs to put limits on the post-shock density and the amount of dust destruction that has taken 
place behind the shock front, as well as put a limit on the dust-to-gas mass ratio in the ambient 
medium, which we found to be a factor of several times lower than the standard value for the LMC 
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of ~ 0.25% (Weingartner & Draine 2001, hereafter WDOl). These limits, however, were based on 
only one IR detection, at 2 4 urn, with uppe r limits placed on the 70 fim detection. Both remnants 



were detected with A^arz (|Seok et al.ll2008|) at 15 and 24 fim (0519 was also detected at 11 /im). 
Seok et al. applied single-temperature dust models to Akari data, deriving warm dust masses of 
8.7 xlQ-^ Mq and 3.6 x 10~^ Mq in 0509 and 0519, respectively. With full spectroscopic data, 
we can place much more stringent constraints on the dust destruction and dust-to-gas mass ratio in 
the ISM. We also explore alternative dust models, such as porous and composite grains. 

Additionally, we examine archival data from the Reflection Grating Spectrometer (RGS) on- 
board the XMM-Newton X-ray observatory. The high spectral resolution of RGS allows us to 
measure the strength of lines in X-ray spectra. We can use post-shock densities derived from IR 
fits to predict the strength of lines arising from shocked ambient medium, and can thus derive 
relative strengths of ejecta contributions to oxygen lines. Doing this requires knowledge of the 
shock compression ratio, r, defined as nn/no (where uh and uq are postshock and preshock hy- 
drogen number densities, respectively), which for a standard strong shock (Mach number ^ 1) is 
4. However, shock dynamics will be modified ( Jones & Ellison|[r991 ) if the shock is efficient at ac- 



celerating cosmic rays, and r can be increased by a factor of several. We use several representative 
values of r in modeling X-ray line strengths. 



2. Observations and Data Reduction 

We mapped both objects using the long wavelength (14-38 /im), low-resolution (AA/A 64- 
128) (LL) spectrometer on Spitzer's IRS (Program ID 40604). For 0509, we stepped across the 
remnant in seven LL slit pointings, stepping perpendicularly 5.1" each time. This step size is half 
the slit width, and is also the size of a pixel on the LL spectrograph. We then shifted the slit po- 
sitions by 56" in the parallel direction and stepped across the remnant again, ensuring redundancy 
for the entire object. This process was repeated for each of the two orders of the spectrograph. 
Each pointing consisted of two 120-second cycles, for a total observing time of 6720 s. 

For 0519, the process was identical in terms of number of positions and step sizes in between 
mappings, but each pointing consisted of four 30-second cycles, for a total observing time of 3360 
s. The difference was due to the fact that 0519 is several times brighter in the wavelength range of 
interest, and we did not want to saturate the detectors. 

To estimate the uncertainties in the data, we fit a second-order polynomial to the line-free 
wavelength regions (~ 21 — 32 /im) in the spectra where the signal-to-noise was adequate for 
modeling. The standard deviations of the data points in these regions were used for fitting of 
dust emission models (see Section |4~TI) . We obtained a values of 0.830 and 2.53 mJy for 0509 and 
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0519, respectively. 

The spectra were processed at the Spitzer Science Center using version 17.1 of the IRS 
pipeline. We ran a clipping algorithm on the spectra which removes both hot and cold pixels that 
are more than 3-a away from the average of the surrounding pixels. To extract the spectra, we used 
SPICE, the Spitzer Custom Extraction tool. Once all spectra were extracted, we stacked spectra 
from the same spatial location, improving the signal-to-noise ratio of the sources. For background 
subtraction, we use the off source slit positions that come when one of the two slit orders is on the 
source. In the end, we have seven (overlapping) background subtracted spectra for each remnant. 
As we show in Section [3l this allows us to do spatially resolved spectroscopy, despite the fact that 
the remnants are only ~ 30" in diameter. 

We processed archival XMM-Newton RGS data from the XMM Science Archive with version 
8.0 of the Science Analysis Subsystem (SAS) software for XMM. RGS provides high-resolution 
(AA ~ 0.04 A) X-Ray spectroscopy in the wavelength range from 5-38 A. B0509-67.5 was ob- 
served on 4 July 2000 (Obs. ID 01 1 1 130201, PI M. Watson) for 36 ksec. Of these 36 ksec, only 
the last ~ 2500 s were affected by flaring. We discarded these and use the remaining 33.5 ksec for 
our analysis. We use the spectra from both RGS 1 and 2, and bin the data using the FTOOL grppha 
to a minimum of 25 counts per bin. Both RGS orders contained ~ 8000 counts after time- filtering. 
B0519-69.0 was observed on 17 September 2001 (Obs. ID 0113000501, PI J. Kaastra) for 48.4 
ksec. After time filtering, we obtained ~ 35 ksec of useful data, with ~ 12000 counts in both 
RGS 1 and 2. We also grouped these by a minimum of 25 counts per bin. Since RGS is a slitless 
spectrometer, spatial information is degraded for extended sources, and the spectrum is smeared by 
the image of the source. In order to model RGS spectra, the response files generated by SAS must 
be convolved with a high-resolution X-ray image. For that purpose, we used archival broadband 
Chandra images from 2001 (Obs. ID 776, PI J. Hughes), along with the FTOOL rgsrmfsmooth, to 
produce new response matrices. 



3. Results 
3.1. 0509 

In the top left panel of Figure [H we show the MIPS 24 yum image of 0509, with overlays as 
described in the caption. Immediately obvious is the large asymmetry in 0509, where the remnant 
brightens from the faint NE hemisphere to the brighter SW. Dividing the remnant in half and 
measuring the 24 jim fluxes from each half, we find a flux ratio of five between these two regions. 
At X-ray wavelengths, we measure a much more modest ratio of 1.5 from archival broadband 
Chandra images. Using IRS, we separated the remnant into 7 overlapping regions from which 
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spectra were extracted. For our analysis, we chose two regions (shown on Figure [T]) that do not 
overlap spatially and provide adequate signal to noise spectra. These regions roughly correspond 
to the bright and faint halves of the remnant, and the spectra are shown in Figure [2l In Figure [3] and 
Figure IH we show both of these spectra separately, with models overlaid (models are discussed 
in Section ID). The ratio of integrated fluxes from 14-35 fim for the two regions is ~ 2.5. This 
differs from the factor of 5 measured from the photometric images because the slits are not exactly 
aligned with the bright and faint halves of the remnant, and because of the differences in the MIPS 
and IRS bandpasses. 

A few obvious features immediately stand out about both spectra. First, there is no line 
emission at all seen in either spectrum. This is not unex pected, since the shocks in 0509 are some 



of the fastest SNR shocks known, at > 5000 km s~^ (|Ghavamian et al.l l2007h . and there is no 
optical evidence for radiative shocks. Second, though both spectra show continuum from warm 
dust, there are obvious differences in the spectra. An inflection around 18 /im can be seen in the 
spectrum from the bright half, while this feature is not as clear in the faint spectrum. We attribute 
this feature to the Si-O-Si bending mode in amorphous silicate dust. In Section |71 we will explore 
reasons for the different spectral shapes and the differences in brightness between the two halves 
of the remnant. 



3.2. 0519 



In the bottom left panel of Figure [U we show the 24 /xm image of 0519, which does not 
show the large scale asymmetries that 0509 does, although three bright knots can be seen in the 
image. These same three knots can be identified in both Ha (|Smith etal.lll99l ') and X-ray images, 
where Chandra broad band data shows that the three knots are prominent in the 0.3-0.7 keV band 
(IKosenko et al.ll2010l) . In Paper I, we measured the total flux from the three knots added together 
and found that collectively they represent only about 20% of the flux. Nevertheless, spectra ex- 
tracted from slits that contain a bright knot do show differences in continuum slope from those 
extracted where no knots are present. We discuss our interpretation of these knots in Section |71 As 
with 0509, there are no lines seen in the IRS data (see Figure [5]), as shock speeds are also quite 
high in this object. 



4. IRS Fits 



We now turn our attention to modeling the emission seen in IRS with numerical models of 
coUisionally heated dust grains. We follow a procedure identical to that followed in our previous 
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work on these objects; see Paper 1. The heating rate of a dust grain immersed in a hot plasma 
depends on both the gas density (collision rate) and gas temperature (energy per collision). In 
general, the electron temperature of a shocked plasma can be determined from fits to X-ray spectra. 
In the case of these objects, plane-shock models without coUisionless heating of electrons at the 
shock front yield the ter nperature. D ust heating models do depend on this temperature, but the 
dependency is not large (lDweklll987|) . Gas density is a much more sensitive parameter for dust 
heating models, and is determined by the temperature of the grains. We regard postshock gas 
density as a free parameter in our model, and fine tune the density to match the observed IR 
spectra. 

We take into account sputtering by ions, which destroys small grains and sputters material 
off large grains. We use a plane-shock model which superimposes regions of increasing shock 
age (or sputtering timescale) Tp = Updt, where Up is post-shock proton density, while keeping 
temperature behind the shock constant. Inputs to the model are an assumed pre-shock grain-size 
distribution, grain type and abundance, proton and electron density, ion and electron temperature. 
The code calculates the heating and sputtering for grains from 1 nm to 1 fim, producing a grain 
temperature and spectrum for each grain size. The spectra are then added in proper proportions 
according to the post-sputtering size distributio n to produce a model s pectrum that can be compared 
with observ ations. We use spu ttering rates from^ Nozawa et aL ( 2006 ) with enhancements for small 
grains from I Jurac et al.l (|1998[) . We assume that sputtering yields are proportional to the amount 
of energy deposited by incoming pa rticles, accounting for partial tra nsparency of grains to protons 
and alpha particles at high energies dSerra Diaz-Cano & JoneslboOS ). particularly relevant for 0509 
and 0519. We assume that grains are compact spheres, but we also report results and implications 
if grains are porous (i.e. contain a nonzero volume filling fraction of vacuum). 

In general, our method is to fix as many of t he input pa r amete rs as possible, based on what 
is known from other observations. For instance, iRest et al.l (|2005b used optical light echoes to 
constrain the ages of both remnants, yielding ages of ~ 400 yrs. for 0509 and ~ 600 yrs. for 0519. 
Shock velocities ar e known approxirnately f or both objects from measurements of ion temperatures 
(|Tuohy et al.lll982l : iGhavamian et al.ll2007|) . In 0509, the standard shoc k jump conditions for a 
7000 km s~^ shock, the assumed speed for the NE half of the remnant (IGhavamian et al.ll2007|) 
with no ion-electro n equilibration would produce proton temperatures of ~ 115 keV. Recently, 
Helder et al.l (|2010|) used VLT spectroscopy to detect a broad Ha component from the rim of 0509, 
determining the shock speed to be ~ 6000 km s^^ in the NE. They give several possibilities for the 
proton temperature in the NE, and note that energy from escaping cosmic-rays limits the post-shock 
proton temperature to a value of < 0.85 of that expected from the strong shock jump conditions. 
Adopting this value of 0.85, along with a shock speed of 6000 km s^^ and no electron heating 
gives a post-shock proton temperature of <70 keV. We use this value as an input to dust heating 
models, except where noted otherwise (see Section 17121) . This upper limit on proton temperature 
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makes any fit to the density for 0509 a lower limit, and we discuss the effects of changing the 
proton temperature in Section |T!2l 



The shock speed in 0519 is a bit more uncertain. iGhavamian et al.1 (|2007l) q uote a range of 
2600-4500 km s^ between the limits of minimal and full ion-electron equilibration. iKosenko et al. 



(|2010t) measured line widths of X-ray iron lines seen in XMM-Newton RGS spectra to obtain 
a velocity of 2927 km s^^. We adopt a value of 3000 km s^^ for this work and use the standard 
shock jump proton temperature assuming no electron-ion equilibration, which is 21 keV. The actual 
proton temperature may differ from this estimate, but the fitting of dust spectra in the IR is only 
weakly dependent on the proton temperature (see Section iT?!]) . 

The overall normalization of the model to the data (combined with the known distance to the 
LMC) provides the mass in dust. Since our models calculate the amount of sputtering that takes 
place in the shock, we can then determine the amount of dust present in the pre-shock undisturbed 
ISM that has been encountered by the blast wave. While we are sensitive only to "warm" dust, 
and could in principle be underestimating the mass in dust if a large percentage of it is too cold 
to radiate at IRS wavelengths, it is difficult to imagine a scenario where any significant amount 
of dust goes unheated by such hot gas behind the high- velocity shocks present in both remnants, 
particularly given the upper limits at 70 fim reported in Borkowski et al.l (2006) (and the fact that 
there are no radiative shocks observed anywhere in these remnants). 



4.1. Model Uncertainties 

The only fitted parameters in the model are the post-shock proton number density, which 
affects the shape of the model spectrum, and the radiating dust mass, which affects the overall 
normalization. All other variable parameters (e.g. electron density, sputtering timescale) depend 
directly on the density. We use an optimization algorithm to minimize the value of in the regions 
of the spectra with high signal-to-noise (for both remnants, ~ 21 — 32 fim) to find the best- fitting 
value of the density. We then calculate upper and lower limits on the fitted density as the 90% 
confidence intervals in space (where confidence intervals correspond to Xmin + 2.71), using the 
standard deviations listed in Section 2. The errors on both the density and dust mass present are of 
order 25% for 0509 and 10% for 0519. Larger errors in 0509 are due to a poorer signal-to-noise 
ratio in the spectrum. We assume that the plasma temperature is fixed at the values listed above, 
but again note that the dependency on proton and electron temperature is small (see Section |T2l). 

There are also several sources of uncertainty in modeling emission from warm dust grains that 
are hard to quantify statistically. First and foremost, we are limited by both the spatial resolution of 
Spitzer for such small SNRs and the low signal-to-noise ratio of the spectra. We use a plane-shock 
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model with constant plasma temperature to model dust emission, which is only an approximation 
to a spherical blast wave. Estimates of swept gas mass assume constant ambient density ahead of 
the shock. Sputtering rates for such high temperature protons and alpha particles are probably only 
accurate to within a factor of 2 (but see Appendix A), and we use the simple approximation that 
the sputtering rate varies according to the amount of energy deposited into a grain of given size 
and composition. Lastly, but possibly most importantly, we do not know the degree of porosity 
and mixing of various grain types. 



4.2. IR Morphologies 

Although we have full spectral mapping of both objects, the 10.5" width of the IRS slit 
severely limits the amount of spatially resolved spectroscopy we can do on remnants that are only 
30" in diameter, as both 0509 and 0519 are. Despite this, we do have several disjoint slit positions 
which isolate spectra from the most prominent features in the images. Because dust radiates as a 
modified blackbody spectrum, emission from a small amount of warmer dust can overwhelm that 
from a larger amount of colder dust, particularly at the wavelengths of interest here. We show 
in Figure [3] and Figure |4] spectra extracted from two regions of 0509, which we label the "faint" 
and "bright" portions of the remnant. Fits to these regions require a density contrast of ~ 4 in 
the post-shock gas, with the higher density required in the "bright" region (higher densities means 
hotter dust, hence more short wavelength emission). Although it might be possible that a density 
gradient of this order actually exists in the ambient medium surrounding 0509 (the angular size of 
30" in the LMC corresponds to a linear diameter of about 7 pc), the implications from this model 
lead to several scenarios which are unlikely and require an appeal to special circumstances. 

In comparing the MIPS image with images of the remnan t at other wavelengths, one can 



clearly see an enhancement in the Chandra broadband image (Warren & HughesI 12004 ) on the 



southwest side of the remnant. This enhancement is mostly present in the energy range containing 
Fe L-shell emission lines. The Ha image, shown in Figure [H shows a uniform periphery around 
the remnant, with the exception of a brightness enhancement in the SW, relatively well confined to 
a small region. To compare the images, we convolved the Ha image to the resolution of the 24 jjLm 
image using the MIPS 24 /zm PSF The result is shown in Figured The images are morphologi- 
cally similar, and from these comparisons we conclude that there is not an overall NE-SW density 
gradient in the ISM. Rather, the ISM is mostly uniform except for the SW, where the remnant is 
running into a localized region of higher density. To obtain the conditions for the whole of the 
remnant, we fit a model to the "faint" region, freezing all parameters to those reported in Table 1 
and allowing post-shock density to v ary (which also causes the spu ttering timescale to vary). If we 



assume a standard LMC dust model flWeingartner & Drainell200l|) . we get a post-shock density of 
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uh = 0.59 cm ^. The reduced value for this model was 1.04. Values in Table [21 correspond to 
this "faint" region. 

We can use IRS to estimate the spatial extent of the bright region. First, by looking at the 2-D 
spectra from the LL slit located directly on top of the brightest region in the MIPS 24 /im image, we 
fit a gaussian to the spatial profile (indicated by an arrow in Figure |7]) of the emission between 16 
and 17 fim, and found a FWHM of the brightness profile of 12.7". This wavelength is the shortest 
at which we could achieve a good signal-to-noise ratio for this region. The Spitzer resolution at 
16 /xm is about 4.75". In the direction perpendicular to the slit positions, we had 7 slit pointings, 
each separated by 5.1". Extracting spectra from each of these (overlapping) positions, subtracting 
the background, and integrating over the IRS bandpass gave us 7 measurements of the brightness 
as a function of position. The locations from which these spectra were extracted are shown as 
diamonds in Figure |71 and the resulting normalized intensity function is shown in Figure [8l The 
brightness peaks at position 3, the position we used to examine 2-D spectra, as described above, 
and falls off relatively quickly on both sides away from that position. 

Our assumption that the bright region of 0509 is the result of a more localized enhancement 
in density implies that the spectrum extracted from that region is then a composite of emission 
from the uniform parts of the remnant and from the small, denser region. Spitzer does not have 
the resolution required to isolate this dense region, but we can make a crude isolation by assuming 
that the two slit positions in Figure [T] cover an equal surface area on the remnant, and subtracting 
the faint spectrum from the bright. The residual spectrum is of less than ideal signal-to-noise. 
Nevertheless, a fit to this spectrum implies a gas density of ~ 6.4 cm~'^, or about an order of 
magnitude higher than the rest of the remnant. A small region of hot dust outshines a more massive 
region of cooler dust at 24 fim, which would explain the factor of 5 ratio in the flux between the 
two sides. X-rays from this object are dominated by ejecta, not swept-up ISM, and there are several 
factors to consider beyond density when considering the flux of Ha coming from a region. Higher 
resolution observations at all wavelengths will shine more light on this issue. 

For 0519 we have a similar scenario, except that instead of one bright region, we have three 
bright knots. As noted previously, these knots correspond spatially with knots seen in both X-rays 
and Ha (see Figure [T]). We adopt an identical strategy here, isolating spectra from slit positions 
that do not overlap with one of the bright knots. As it turns out, there is only one slit position in this 
object that is nearly completely free of emission from a knot, a slit position that goes directly across 
the middle of the remnant. The spectrum from this region is shown in Figure [51 We assume that 
the conditions within this slit position are indicative of the remnant as a whole. Using parameters 
found in Table 1, we obtain a post-shock density of uh = 6.2 cm^^ (reduced = 1-13). We do not 
have adequate signal-to-noise to isolate the bright knots individually. However, we approximate 
the density required by noting that a fit to the spectrum extracted from the entire slit position that 
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contains the brightest of the three knots (the northernmost knot) requires a density that is higher 
by a factor of ~ 2. Given that this spectrum contains emission from both the uniform parts of the 
remnant and the bright knot, it is likely that the density in the knot itself is perhaps a factor of ~ 3 
higher that the more uniform parts of the remnant, which is roughly consistent with our results 
from Paper I. 



5. Dust Grain Porosity 

While it is typically assumed in dust grain models that gr ains are solid bodi es composed of 



various materials, this is likely not a physically valid model (see lShen et alj (120081) . and references 
therein). It is possible to relax the assumption that grains are compact spheres with a filling fraction 
of unity. We have updated dust models to include the porosity of dust grains (where porosity, V, is 
the fraction of the grain volume that is composed of vacuum, instead of grain material), as well as 
allowing for grains to be composite in nature, i.e. a grain made u p of silicates, grap hite, amorphous 



carbon, and vacuum. The details of this model are discussed in iWilliamsl (120101) . We summarize 
here the important results. Perhaps the most significant effect of introducing porosity to grains is 
that the grain volume per unit mass increases, since the overall density of the grain is lowered. 
As a result, grains are more efficient radiators of energy, and less dust mass in total is required to 
reproduce an observed IR luminosity. 

For a grain of the same mass, the grain is heated to a higher temperature if it is solid than 
if it is porous. As a result of the equilibrium temperature being lower for porous grains, more 
heating is necessary to fit the same spectrum with a porous grain model than a solid grain model. 
For a fixed plasma temperature, this can only be achieved by increasing the post-shock density. In 
summary, increasing the porosity of grains has the dual effects of raising the necessary postshock 
density required to heat grains to an observed temperature, while lowering the total amount of dust 
required to reproduce an observed luminosity. 

What remains is to choose a grain "recipe" to use to fit the data. Since there are few constraints 
from the ISM at this time, we have adopted here two values of V, 25 and 50%, to explore the effects 
of grain porosity of fitting IR data. Our ''compact" model, which we used to derive results given 
above, assumes separate populations of solid silicate and graphite grains with ap propriate size 
distrib utions for the LMC (WDOl). Our ''porous" grain models are based on those of IClayton et al. 



(120031) . For the 25% porous model, we assume that grains are composed of 25% vacuum, 50% 
silicate, and 25% amorphous carbon. The 50% porous grains contain 50% vacuum, 33.5% silicate, 
and 16.5% amorphous carbon, distributed from 0.0025 to 1.5 /im. The narrow bandwidth of the 
IRS data does not allow us to determine which model fits the data better, as we obtained equally 
good fits from all models (albeit with a different value of density), and report results from all. A 
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caveat to using these porous grain distributions is that both of them were optimized for the Galaxy, 
and not the LMC. However, to the best of our knowledge, a size distribution appropriate for porous 
grains in the LMC has not yet been developed. 

The results are as follows: For 0509, the best-fit V = 25% model had a post-shock density, 
uh, of 1.1 cm~^ (x^ = 1.05), while the V = 50% model required an uh of 2.2 cm~^ (x^ = 1.07). 
For 0519, the best-fit densities were 9.5 (x^ = 1.08) and 17 (x^ = 1.07) cm'^ for P = 25 and 50%, 
respectively. In all cases, we assume LMC abundances, such that rig = I.Itih- 



6. X-ray Modeling of RGS Data 

Although dust models in the IR are a powerful diagnostic of the post-shock gas, they are 
insensitive to pre-shock gas conditions, which determine the total swept mass. Unlike many SNRs, 
both of these objects have relatively well-defined ages and distances. The post-shock densities 
are related to the pre-shock values by r, which is 4 in the case of an unmodified, strong shock, 
but may be greater than 4 in the case of cosmic-ray modification of the blast wave. For a given 
r, we can calculate the pre-shock density, and thus the total amount of swept up material. The 
product of post-shock electron density and total gas mass overrun by the blast wave is proportional 
to the emission measure (EM) of the X-ray emitting plasma. When combined with the ionization 
timescale (Tj = n^dt) and electron temperature of the gas, we have all the information necessary 
to generate a model X-ray spectrum for the shocked ambient medium, which we can compare 
directly to observed X-ray spectra. 

The difficulty in this approach lies in disentangling X-ray emission from the shocked ambient 
medium from that arising from reverse- shocked ejecta. For large remnants, this could be done 
spatially using CCD spectra, but both 0509 and 0519 are ~ 30" in diameter with the ejecta not 
well sep arated from the shocked ISM, rendering such an approach impossible. For 0509, previous 



studies (jWarren & Hughesll2004l : iBadenes et al.ll2008|) have used Chandra CCD spectra to model 



emission from the remnant, but these studies are subject to the inherent difficulty in using poor 
spectral resolution CCD spectra to disentangle ISM and ejecta emission. An alternative approach 
is to take high-resolution spectra from grating spectrometers, such as RGS on XMM-Newton. If 
lines can be identified as arising only from shocked ambient medium, then model fits to these 
lines could be directly compared with predicted line strengths from the model described above. 
This can be most eas ily done with N lines, as virtually no N is produced in Type la explosions. 



Kosenko et al.l (l2008b fit the 24.8 A N Ly-a line in the spectrum of 0509, but we find that this line 
is very weak and unconstraining to the fitting procedure used here. Here we generate model O 
Lya and Ka lines visible in RGS spectra and compare them directly to observed line fluxes. In 
Figure m we show RGS spectra from both remnants. 
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6.1. 0509 



As seen in Figure HI Fe L-shell lines dominate the spectrum of 0509, but strong lines from 
both H and He-like O are clearly visible in both cases at 18.97 and 21.7 A. In order to measure 
the strength of observed O lines, we fit the data with a two-component model. One model, a 
vpshock model containing only iron, was fit to the data between 14-18 A, which are dominated 
by iron lines. We then added a second vpshock model at fixed normal LMC abundances (wilms 
abundance model in XSPEC dWilms et al.ll2000|) . with abundances of heavy elements set to 0.4; 
except that N = 0.1) with ionization timescale and normalization free to fit the data between 9-14 
and 18-23 A, regions where the signal-to-noise was good enough to constrain the fit. We use two 
absorption comp onents, one resulting from G alactic absorption with column density, N^, equal to 



5 X 10^° cm"^ ( Staveley-Smith et all 120031). and one from the LM C (at fixed LMC abundances; 



same as above), equal to 2 x 10 cm (iWarren & Hughesll2004|) . These models were used only 



to determine the strength and width of oxygen lines. 

We freeze the electron temperature for the second vpshock model to 2 keV. This is roughly 
the number one gets from models of Coulomb heating of electrons behind a shock of 7000 km s^ 
with a preshock density of 0.25 cm~^. Although we do not know the preshock density a priori, we 
use this as a representative value and note that the dependence of the emission measure of the gas 
on the electron temperature is small. We introduce an additional line smoothing parameter to both 
models to allow for the fact that lines are broadened by Doppler motions. For this model, we find 
that oxygen line fluxes measured from RGS spectra in 0509 are as follows: O Lya = 4.41 x 10^^'^ 



ergs cm ^ s 



The oxygen line widths obtained from 



ergs cm ^ s ^ and O Ka = 5.05 xlO 
fitting are vfwhm = 12350 (11220, 13630) km s^^ 
of Fe lines in the RGS spectrum of this object, obtaining line widths of 1 1500 (10500, 12500) km 



Kosenko et al.l ( 20081) measured the widths 



6.2. 0519 



Our technique for fitting RGS data from 05 19 was identical, except that we extended the fitting 
range out to 27 A to account for the fact that we had better signal-to-noise at long wavelengths in 
this case. We again used tw o absorption components; one for the galaxy with Nm = 6 x 10^° cm~^ 
Jstavelev-Smith et al.ll2003l) and one for the LMC with iV^ = 1 .6 x 10^^ cm ^ faughes et all 1995 



Kosenko et al 



2010l) . Measured line strengths for this remnant are 5.81 X 10 ergs cm s and 



3.91 X 10 ergs cm ^ s ^ for O Lya and O Ka, respectively. We obtain a line wi dth for oxygen 
lines o f vfwhm = 3470 (3170, 4000) km s"^. Using a three-component NET model. Kosenko et al 



(120 10|) measure Fe line widths in this object to be 4400 (4280, 4520) km s ^. The also report that 
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neither Ne nor Mg is required to fit the XMM spectra. We find that normal abundances of Ne and 
Mg from the shocked ISM are allowed by the data, and it is not necessary to exclude them from 
the fitting procedure. 

Oxygen line widths in 0519 are comparable to the assumed shock velocity, while those in 
0509 are significantly higher. In Section l73l we discuss hydrodynamic modeling of 0519, and find 
that given the age of 05 19, a lower shock speed of 2450 km s^^ is preferred to match the observed 
radius of the forward shock. Line widths from a shock of this speed would fall within errors of the 
observed line widths above. 



7. Discussion 

IR spectral fitting gives us a postshock density, which we can then divide by r of the shock 
to get the preshock density. For a standard, "strong" shock, r is pi/po = 4, but modification of 
shocks by cosmic ray acceleration will increase this. An additional complication for the cases of 
these two objects is that they are both only 30" in diameter, which is quite small for the telescope 
resolution at Spitzer wavelengths, meaning that the spectra shown in Figures 2-4 are the average 
spectra over the entire postshock region. To account for this, we computed hydrodynamic models 
(see Section 17.31) and calculated the mean, EM-weighted r over the post-shock region from the 
leading edge of the blast wave to the contact discontinuity. For the case of a standard shock, this 
mean r is 2.65. This is lower than the immediate postshock value of 4 due to the sphericity of the 
blast wave. The postshock density drops behind a spherical shock. 

Taking the example of 0509, dust model fits to IR spectra give a postshock density of n^^^ = 
0.59 cm^^. Dividing this by 2.65, we get an ambient pre-shock density, no, of 0.22 cm~^. The 
known distance to the LMC of 50 kpc, combined with the measured angular size of the remnants, 
gives a shocked volume of 5.7 x 10^'' cm^. Multiplying the pre-shock density by the volume gives 
the value reported in table 2 for the swept-gas mass, 1.47 Mq. Dividing the dust mass inferred 
from fits to IR data, 1 .9 x 10~^ Mq, by this number, we arrive at an ambient dust-to-gas mass ratio 
of 1.3 X 10 the number we report in Table 1 . This dust-to-gas ratio is not folded into the models 
at any point, nor is it ever assumed. We repeat these steps for various values of r, determining 
different swept gas masses and dust-to-gas mass ratios. However, note that the post-shock density 
and total dust mass remain constant for each case, since those were determined from fits to IR 
data. We include effects of sputtering of dust grains by ions in deriving the dust masses above. 
In Figures [TOl and [TTl we show the best fit values of post-shock density and dust mass for various 
grain models. 

X-rays in 0509 and 0519 are coming mostly from reverse-shocked ejecta. The shocked ambi- 
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ent medium may contribute substantially, particularly to continua and O lines, though it is difficult 
to evaluate this contribution based solely on X-ray data. However, because fits to IR data provide 
us with the post-shock density, and the ages and sizes of these remnants are well-known, we can 
use X-ray emission codes to predict line strengths coming from the shocked ambient medium, 
assuming typical LMC abundances. 

In Tables 2 & 3, we show predicted oxygen line strengths as a function of r and grain porosity 
for both remnants. These line strengths were calculated in XSPEC by taking the plane- shock mod- 
els as described in the section above and varying the ionization timescale and emission measure 
appropriately for each case. Real ages of 400 and 600 years for 0509 and 0519, respectively, were 
divided by 3 in calculating Tj as a c orrection factor for applying a plane-shock model to a spherical 



blast wave (iBorkowski et al.ll2001|) . EM was calculated as rieTioV, where rie is post-shock electron 
density, uq is pre-shock density (= where V is total volume enclosed by the blast wave (ap- 
prox. 5.7 X 10^^ cm^ for both remnants). We chose values for r of 4, 8, and 12, which correspond 
to mean EM-weighted r values of 2.65, 4.8, and 6.8. 

In the Tables, we quote our predicted O line strengths from the blast wave as a fraction of 
the observed strengths from RGS spectra. Values < 1 imply a contribution to O lines from ejecta, 
while values > 1 are inconsistent with observations and are ruled out. For 0519, this rules out all 
models of porous grains. However, this requirement also rules out compact grain models that are 
unmodified from "standard" shock jump conditions. Only those shocks with significantly higher r 
are possible based on X-ray line strengths observed. 

For 0509, only one case is ruled out, that of highly porous {V = 50%) grains heated behind 
a shock with standard jump conditions. If grains are compact, then virtually all O line emission 
seen in RGS spectra wou ld have to come from ejecta. Recent near-IR observations of Type la 



SNe (jMarion et al.l 120091) suggest that the entire progenitor is burned in the explosion and that 



O and Ne are byproducts of carbon burning, meaning that strong X-ray lines from oxygen in 



the ejecta may be unlikely. The favored explosion model for 0509 from iBadenes et al.l (|2008|) 
contains 0.04 Mq of unburned oxygen, the least of any of the explosion models considered for 
this remnant. Nevertheless, even this relatively small amount of ejecta oxygen could be enough 
to produce observed O lines. Model line ratios require at least some contribution from reverse- 
shocked ejecta to account for observed O hya strength. 



Helder et al.l(l2010t) find that at least some of the energy of the shock in 0509 must be deposited 
into cosmic rays. Of course, if cosmic -ray pressure is indeed non-negligible, then the standard 
case no longer applies, and cosmic-ray modified shocks would be required. Our assumed proton 
temperature of 70 keV already implies modification from the standard shock jump conditions, 
which would predict kTp ~ 115 keV. However, the dependency of the dust heating models on 
this number is small (see Section 17121) . The relationship between cosmic-ray pressure and shock 
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compression ratio is model dependent (|Vink et al.ll2010|) and beyond the scope of this paper, but 
models in Table 1 with r > 4 would be preferred. 

The data do not allow us to firmly distinguish between compact and porous grain fits, as both 
can be consistent with the data under the right circumstances. Instead, we summarize implications 
for each model. If dust grains in the ISM are compact, the oxygen lines visible in X-ray spectra 
of 0509 are dominated by contributions from the ejecta, for all values of r. The inferred preshock 
density is inversely proportional to r (see Table [21). A "standard" r value of 4 yields a preshock 
density of nn = . 22 cm ~'^, significantly higher th an either the value of 0.05 cm^'' quoted by 
Warren & HughesI hoo4) or that of 0.02 cm^ froml Tuohv et all (Il982h . but below the 0.43 cm^ 



of Badenes et al. 



(|20of In 0519, it is possible for nearly all of O lines to be coming from shocked 
ambient medium; however, the minimum r value required to match the O Ka line strength observed 
is 12. For this value of r, a small ejecta contribution to the O hya line is still necessary for line 
ratios to match. 

If grains are porous, a lower ejecta contribution to oxygen lines observed in 0509 would be 
permitted, even for values of r that are closer to the standard shock value. However, in 0519, much 
higher r values would be required to bring model oxygen line fluxes down to observed levels, 
even with no contribu tion from the ejecta. We require significantly less oxygen in the ejecta than 
Kosenko et al.l (|2010|) . who report an O ejecta mass of 0.2-0.3 Mq. Minimum values for r would 
be 20 and 34 for V = 25% and 50%, respectively. 



7.1. Radio Limits on Compression Ratio 

We can set an upper limit to r using the observed diffuse synchrotron brightness of the LMC 
near 0509 and 05 19, and requiring that the brightening by compression of the ambient synchrotron 
emissivity not exceed the observed radio brightnesses of 0509 and 0519. Compression will boost 
the energy of electrons and also increase the mean magnetic field strength in the emissivity ju- If 
the electron energy distribution is N{E) = KE^^ between energies Ei and E^, then the emissivity 
is given by 

= Cj{a)KB^+''v~'' (1) 

where a = (s — l)/2 and Cj(0.5) =4.85x10"^"^ cgs. (In general, Cj = C5 (a) (2ci)" in the notation 
ofPacholczyk [1970], with ci = 1.82 x lO^^ and 05(0.5) = 1.37 x IQ-^s.) Taking a = 0.5 (s = 2) 
for simplicity, we find that the electron energy density Ue obeys 

Me = y ' KE-'EdE = K\n (^^^ oc Pe oc p^J^ oc r^l^. (2) 
In the absence of turbulent effects and magnetic-field amplification,the tangential component 
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of B will be increased by r while the radial component will remain the same, raising the magnitude 
of i? by a factor < r. For a spherical remnant encountering a uniform magnetic field, it is 
straightforward to show that the mean amplification factor (r^) is 0.79r for r > 6. Thus, after 
compression, the synchrotron emissivity interior to the remnant is given by 

= r4/3(o.79r)i-5 = OJOr^^/^ 

For a mean diffuse surface brightness upstream of Sq W m~^ sr^^, the mean emissivity in the 
LMC just depends on the line-of-sight depth of synchrotron-emitting material in the LMC. Calling 
that quantity /, we have (j;/(up)) = Sq//. Compressing this as above predicts the post-shock 
surface brightness for a SNR of radius R, /snr^ to be about 

W = J.(in)i? = 0.70r''/'Mont)R = 0.70r''/' (^^^ R (4) 

/ T \ 6/17 / , N 6/17 

^,„.,13(^) (1) . 

(6) 

The diffuse radio brightness of the LMC has been mapped with ATCA and Parkes (Hughes et 
al. 2007). From the image at 20 cm (resolution 40"), we estimate surface brightnesses near 0509 
and 0519 of about 1 and 2 mJy/beam respectively, or about 1 x 10^^^ and 2 x 10^^^ W m^^ sr^^ at 1 
GHz, respectively (using a mean spectral index of the diffuse emission of 0.3; Hughes et al. 2007). 
We take / ~ 1 kpc. 



For 0519, a high-resolution image is available (|Dickel & Milnelll995 l). which shows a shell 
peak brightness of about 6 mJy/beam at 20 cm, giving (for a = 0.5 and a FWHM of 4.8") Jsnr ~ 
1 X 10^^^ W m^^ sr^\ and rmax ~ 26. There is no high-resolution image of 0509 available; using 
the mean remnant surface brightness (~ 5 x 10^^° W m^^ sr^^ from data in Mills et al. 1984) 
slightly overestimates rmax, but gives a value of 25. Larger values of r would give higher radio 
brightnesses than we observe. 



7.2. Effect of Decreased Ion Temperature on IR Fitting 



If shocks are indeed modified by cosmic-ray accleration, both r and the pla sma temperature 
behin d t he shock wou l d be a ltered from the standard strong shock jump conditions (|jones & Ellison 
199 iHelder et al.l (|2010l) report evidence in 0509 for decreased proton temperatures due to 
cosmic-ray modification of the shock. In this work, we find that higher than normal values of 
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r are also necessary to avoid overpredicting oxygen lines in the X-ray spectra of 0519; such a 
modification would likewise lower proton temperatures from the value assumed. Dust model fits 
to IR data are dependent only on these post-shock conditions, and are insensitive to the degree of 
modification at the shock. A detailed examination of this relationship is beyond the scope of this 
paper, but we can make some general comments about the effect that lowering the proton temper- 
ature behind the shock would have on fitting the data. In the case of most SNRs, heating of grains 
by protons can be neglected, as heating is dominated by collisions wit h electrons. Ho wever, for 
remnants with very fast shocks, proton heating becomes non-negligible dWilliamslllOlOh . 

In Table[2l we use a proton temperature of 70 keV for modeling 0509. We take the case of 25% 
porous grains as an example (where uh was fit to be 1 . 1 cm^^). Raising the proton temperature by a 
factor of ~ 50% to 115 keV, the temperature expected from a 7000 km s~^ shock with no electron- 
ion equilibration lowers the fitted density to 1.0 cm an effect of only about 10%. Lowering by 
a factor of ~ 10 to Tp = 7.5 keV raises uh to 1.3 cm^'^. Finally, setting the p roton temperature 



to Tp = Te = 2 keV (a lower limit; although unphysical in this object given that Gha vamian et al. 



(|2007h report detection of broad Ly/3 UV line with width 3710 km s ^) requires uh = 2.0 cm ^. 
These effects are non-negligible, but are not enough to affect the overall conclusions stated above. 



7.3. Hydrodynamic Modeling of 0519-69.0 

We can derive a conservat ive lower limit on t he pre-shock density in 0519 using one of two 
methods. The Ha fluxes from iTuohy et al.l (|1982h give a firm lower limit on tiq of 0.2 cm at 



our assumed shock speed into a fully neutral medium, and assuming 0.2 Ha photons per atom 
entering the shock. The radio upper limit to r (see Section 1711) of 25, which translates to a mean 
r of ~ 12.5, combined with the minimum post-shock density of 6.2 cm^'^ yields a lower limit 
on no of 0.5 cm~^. An upper limit on the pre-shock density in 0519 is given in Table 3 as 0.91 
cm~^, which i s a factor of ~ 2.5 lower than the pre-shock density of uq = 2.4 cm~^ quoted by 
Kosenko et al. IfiolO), who derive their ambient density by fitting a six-component NEI model to 



the X-ray spectra from Chandra and XMM-Newton. In their model, most of the oxygen comes 
from the ejecta, which requires a very high ionization timescale of 8 x 10^^ cm^'^ s (a factor of 
~ 20 higher than our value reported in Table 3) and thus a very high EM (2.36 x 10^^ cm^) for the 
shocked CSM component, implying a higher pre-shock density. 

There is roughly a factor of 2-3 between our lower and upper limits to the pre-shock den- 
sity. To explore both of these limits, we use the Virginia Hydrodynamics (VH-1) time-dependent 
hydrodynamic numerical code in one dimension. VH-1 is a hydrodynamics code which solves 
the Eulerian equations of fl uid mechanics using the Piecewise Parabolic Method, described in 



Colella & Woodward! (|1984h . We assume spherical symmetry on a numerical grid with 512 zones. 
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We modify r for the forward shock by varying the adia batic index, 7 (= | for a normal shock), of 



the shocked ambient medium ( Blondin & EUisonlboOl ). Softening the equation of state effectively 
makes the sho cked plasma more compressible , yielding r > 4. We use the exponential ejecta den- 
sity profiles of Dwarkadas & Chevalier ( 1998 ^ and assume an explosion energy of 10^^ ergs with 
an ejected mass of 1.4 Mq. In Figure [T2l we plot the evolution of the forward and reverse shocks 
with time for varying values of 7 ^ corresponding to immediate p ostshock r values of 4, 8, and 12, 
using identical scalings to that of IDwarkadas & Chevalied (|l998h (compare with Figure 2a of their 
paper). 

With the mean post-shock density fixed at 6.2 cm~'^ (the gas density required to fit the spec- 
trum assuming compact grains), we examined r values of 12 and 25, corresponding to 7 values of 
1.18 and 1.085, mean r in the post-shock region of 6.8 and 12.5, and pre-shock densities of 0.91 
and 0.50, respectively. 

The lower density case, with r = 25, reaches a radius of 3.6 pc in 550 years with a shock 
velocity at that time of 2950 km s^^ While this age is roughly consistent with the known age of 
this object from light echoes, the shock velocity is too high to be consistent with the measured 
width of oxygen lines. This shock speed would imply a width, vfwhm, of 5660 km s^^. The 
measured value was ~ 3500 km s~^. The higher density case requires 630 years to reach the 
observed radius, with a current shock speed of 2450 km s^^ The line width implied from this 
shock speed is 4490 km s^^ While this is still higher than our observed line width, it is much 
closer to being within the quoted errors, and is clearly favored over the lower density case. 

Errors in determining the post-shock density from IR fits are about 20%, but this has almost 
no effect on the inferred pre-shock density allowed by not overpredicting oxygen lines in RGS 
spectra. However, this number is sensitive to the soft X-ray absorbing column density used in the 
models. If the LMC absorption component, detailed in Section [6^ is increased, a higher pre-shock 
density will be allowed before oxygen lines are overproduced. To be within the errors of the line 
widths measured for oxygen lines, the LMC A^^ v alue would only need to be modestly increased 
from 1 .6 to 1 .7 X 10^^ cm~^. iKosenko et al.l (|2010l) report an A^^^ of 2.6 x 10^^ cm~^, which would 
easily accommodate a higher density. Uncertainties in the exact distance to the remnant are also 
non-negligible here. We assume a distan ce to the LMC of 50 kp c, but recent measurements of the 
distance to the LMC range from 47.8 kpc (iGrocholski et al. II2OO7I) to 5 1 . 1 kpc iKoerwei]l2009l) . The 
lower of these two would increase the observed radius of the remnant from 3.6 to 3.77 pc, which 
would bring the observed line width within errors of the shock velocity from hydro modeling. The 
unknown position of 0519 within the depth of the LMC may also be important. 
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7.4. Dust-to-Gas Mass Ratio 



Of particular note in all results presented in Tables 2 and 3 is that given the amount of gas 
that these remnants have swept up at this stage of their evolution, the amount of dust observed 
by Spitzer is lower than what is expected, even accounting for the effects of sputtering (although 
there are a few values in Table 2 that are acceptable; namely, those with a high r). The amount 
of dust sputtered is ~ 10% for compact g rains in 0509, and 34% for compact grains in 0519. 
UV studies of grain absorption in the ISM (|Weingartner & Drainell200l|) yield a dust-to-gas mass 
ratio in the LMC of 2.5 xlO~ ^. Recent wide-field far-IR observations of the LMC give a value 
of ~ 3.5 X 10 for this ratio ( Meixner et al.l 2010). This shortfall is consistent with our previous 
studie s of both LMC and galactic SNRs ( Borkowski et al. 20061: ' ' 



Wi 



2007t) . Similar results have been found for SN 1987A feouchet et al.ll2006 l) in the LMC, although 



iams et al.ll2006l: iBlair et al 



the shocks there probe circumstellar material, and have not ye t reached the an ibient ISM (but see 
Dwek et al. 2008,2010), and for G292.0-I-1.8 Lee et al.ll2009l) and Puppis A dArendt et allboiol) 
in the Galaxy. 

Of all parameters in the model, dust mass is the least model-dependent. As a check of this, we 
examine the case of 0519. Assuming compact grains, we derive a total radiating dust mass of 1.5 
X 1O~^M0 based on model fits to IRS data. We then made fits to the data using a greatly simplified 
model, that of a single grain at a single temperature, ignoring heating, sputtering, and grain-size 
distribution. This model yielded a total dust m ass of a bout 1 .3 x lO'^M©. As a further check, we 



used an approximate analytic expression from lDwekl (|1987l) for the mass of dust, dependent only 
on the temperature of the grain and the observed IR luminosity, obtaining a value of 1.7 x 1O^^M0. 
These three different approaches of decreasing complexity only differ by at most a factor of 30%. 
Taking our minimum value of swept-up gas mass (5.3 Mq), we derive dust-to-gas mass ratios of 
2.4-3.2 xlO""^, an order of magnitude lower than the expected value. Dust-to-gas mass ratios in 
0509 are closer to expected values, but are, in most cases in Table |2l still too low. It may be 
that SNRs that are in regions of low background sufficient to allow a high enough signal-to-noise 
ratio for a detailed study in the mid- and far-IR have such a low background because they are in 
atypical regions of the ISM, creating an observational bias. In addition, dust-to-gas mass ratio 
determinations from SNRs are spot measurements, and are not averaged over long lines of sight or 
large fields of view. 

Porous dust grains require a higher density relative to compact grains to reproduce the same 
spectrum, because they are more efficient radiators of energy than compact grains. For this reason, 
the inferred dust mass from the porous grain models, seen in Table 2, is lower by a factor of ~ 2—4. 
However, for a given r, the inferred gas mass is higher for porous grains, so that the dust-to-gas 
mass ratios for porous models are lower by a factor of ~ 4 for P = 25% and ~ 18 for V = 50%, 
when compared to the compact grain case. It should be noted though, that typical dust-to-gas mass 
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ratios calculated for the ISM in general (such as the 0.25% given above from WDOl) inherently 
assume compact grains. 



One possibility that has been suggested for SNRs of both Type la and CC origin (IGomez et al 



20091 : bunne etaPboOSI) is that large amounts of dust are indeed present, but are too cold (1^ ~ 16 
K) to be detectable at Spitzer wavelengths and are only visible in the sub-millimeter regime. If 
this is true, then the total dust-to-gas mass ratio for the ISM could be quite different from what 
we re port here, which is based only on "warm" dust. This claim has been disputed for Kepler's 
SNR (|Blair et al.ll2007|) based on a detection at 70 /xm and an upper limit at 160 /xm, but we lack 
detection at either wavelength here, and upper limits are unconstraining. Herschel observations 
would provide the key wavelength coverage between the mid-IR and sub-mm necessary to settle 
this issue. A recent Herschel observation of Cas A at wavelengths up to 500 /xm revealed only 
0.075 Mq of dust there, much of which may be newly-formed ejecta dust (IBarlow et al.ll2010|) . 



How cold would dust have to be to escape detection at 70 jim, assuming a standard dust-to-gas 
mass ratio for the LMC? Obviously, the answer depends on the choice of models from Tables 2 & 
3. As an example, we consider the V = model for 0519 with the highest value for r (the only 
model for this remnant that seems to not be ruled out by X-ray lines). The swept gas mass for this 
case is 6.1 Mq, and assuming a dust-to-gas ratio of 0.35% gi ves a dust mass of 2. 1 xlO^^ Mq. 
In order to avoid violating the 70 /xm upper limit of 121 mJy (IBorkowski et al.ll2006|) . post-shock 
dust in this remnant would need to be < 44 K. If we choose the standard r case, the implied dust 
mass would be 5.5 xlO^^ Mq, and the required post-shock dust temperature would be < 36 K. 
Both of these numbers are significantly lower than those fit to IRS data or those predicted by grain 
heating models. Such low temperatures could only be acheived by having the dust reside in very 
dense clumps. Shocks driven into such clumps should quickly become radiative, but no [O III] or 
[S II] emission is seen from either remnant. 



8. Conclusions 

We present mid-IR spectral observations of two young SNRs in the LMC, 0509-67.5 and 
0519-69.0, as well as analysis of archival high-resolution X-ray spectra. By fitting dust heating and 
sputtering models to IR spectra, we can determine post-shock gas density. We derive post-shock 
densities of 0.59 and 6.2 cm^ in 0509 and 05 19, respectively, using a compact grain model. Porous 
grains require a density higher by a factor of a few. By assuming a value for the shock compression 
ratio, r, we can infer the pre-shock density, swept gas mass, strength of oxygen Lya and Ka X-ray 
lines, and dust-to-gas mass ratio for the ambient ISM. Derived values for the pre-shock density 
of the ISM vary greatly depending on the model used and compression ratio assumed, but for the 
sake of comparison, values listed in Tables 2 & 3 are comparable to densities inferred from various 
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line-of-sight HI column densities through the LMC. These column densities vary from 10^° cm~^ 
to a few xlO^^ cm~^, which, assuming an LMC depth of 1 kpc, yields densities of ~ 0.05 — 1.5. 
For standard strong shocks, r is 4, but modification by cosmic rays would increase this by an 
unknown amount. We report values of the quantities above for several plausible values of r. In 
principle, this method could be used in reverse to determine r. This would require X-ray emission 
from shocked ambient medium that is well-separated from ejecta emission; thus. X-ray lines could 
be modeled as solely arising from material at cosmic abundances. We believe this method can be 
used for older, larger remnants in the LMC, such as DEM L71. 

Based on our analysis, a significant fraction of the O line emission seen in the X-ray spectrum 
of 0509 arises from the ejecta. If grains are highly porous, then the ejecta contribution is less, but 
a contribution is still required because the line flux ratios observed still do not match the model 
ratios. Both standard shocks and cosmic-ray modified shocks with higher r values can provide 
acceptable fits for 0509. The data for 0509 are inconclusive in favoring either compact or porous 
grains. Compact grains require nearly all of the oxygen seen in X-ray spectra to be coming from 
shocked ejecta, while porous grains would significantly lower the ejecta contribution, implying a 
higher contribution from the forward-shocked material. 

In 05 19, there is significant evidence for higher than standard compression, and little contribu- 
tion from ejecta to O X-ray lines. A shock speed of < 2300 km s^ is favored for this object from 
hydrodynamical simulations, in order to reproduce the observed size, age, and X-ray line widths. 
Compact grains are favored for this remnant. The standard value of r = 4 is possible only if ab- 
sorption to the remnant is significantly higher than has been reported. The column density from 
the LMC would n eed to be at least 4 x 10^^ cm^^ to bring oxygen lines in the model down to below 



observed values. iGhavamian et al.l (|2007|) derive an upper limit for 05 19 of E(B — V) < 0.11, 



implying an upper limit to the HI column density of ~ 1.6 x 10^^ cm ^ (|Cox et al.ll2006l) 



We derive dust-to-gas mass ratios that are lower by a factor of several than what is generally 
expected in the ISM. Since the general properties of dust in the ISM come from optical/UV absorp- 
tion line studies averaged over long lines of sight, it is possible that there exist large local variations 
on smaller scales. The ratios presented here probe parsec scales. Another potential explanation is 
that the sputtering rate for dust grains is significantly understated in the literature. However, since 
the deficit of dust is about an order of magnitude, sputtering rates would also have to be increased 
by this amount to account for the discrepancy. Also, because sputtering rates typically assume a 
compact grain, further work in this field will be needed to determine if these rates are appropriate 
for porous grains. Extremely high values for r (> 50) could bring dust-to-gas mass ratios up to 
more typical values, but such values are ruled out by radio observations. A possible observational 
bias could exist, in that SNRs that are easy to study in the IR (i.e., well separated from sources of 
IR confusion) could be found preferentially in low-dust regions of the ISM. 
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Studies such as this will benefit greatly from the increases in both spatial resolution and sen- 
sitivity of future generations of telescopes, such as the James Webb Space Telescope. Being able to 
spatially separate the dust spectra right behind the shock from that further inside the shell will be 
crucial to reducing some of the uncertainties listed above. Longer wavelength observations with 
the Herschel Space Observatory would be useful for detecting or placing upper limits on any cold 
dust that may be present in the remnant. 

We thank John Blondin for providing the VH-1 hydrodynamics code and for useful discus- 
sions on interpretation, and the anonymous referee for many useful comments which improved the 
paper. We acknowledge support from Spitzer Guest Observer Grant RSA 1328682. 

Facilities: XMM, Spitzer 
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Table 1. Model Input Parameters 



Object 


Vs (km s-i) 


Te (keV) 


Tp (keV) 


Age (yrs.) 


Ref. 


B0509-67.5 


6000 


2.0 


70 


400 


1,2 


BOS 19-69.0 


3000 


1.5 


21 


600 


1,2 



Note. — 



References: (1) Ghavamian et al 2007, (2) Rest et al. 2005 



Table 2. Predicted ISM Oxygen Line Strengths for SNR B0509-67.5 



V 




no*" 






OLya 


OKa 








Md/Mg 









^■-'^0.46 
^••'^0.46 


0.220;?? 


2.97^1 xlO^ 
2.973;6[ xlO^ 
2.97111 xlO^ 


0-oKill 
o.o36g:8if 


02fi0 0''4 
0.0170:0?i 

0.0140:°^ 


A 0^40.089 

0.04l0:0i 


1 471-82 
1-^' 1.15 
A qO.99 
^■^0.62 

0.68:11 


1.90?i^ xlO- 

1.90?:!^ XlO- 
1.90^:11 XlO- 


-3 
-3 
-3 


1.3^:1 xio- 

2Af l xlO- 

3.2f:i xio- 


-3 

-3 
-3 


25% 


^•^0.9 


o.42g:i 


5.531:^2 xlO^ 


0.32g:il 


n 190.21 


0.390:61 


2.80i:i 


9.37i2s2 xlO- 


-4 


3.3f:i xio- 


-4 


50% 


9 92.7 








o.62g:i 


^•^^0.94 


5 5^6.78 


4.iyi XlO- 


-4 


lAl^j xlO- 


-5 



"^Hydrogen number density in cm"^. 
''Ionization timescale in cm"^ s. 
'^Emission measure (EM) in 10^^ cm-^. 
''Swept gas mass in Mq. 
•^Dust mass from fit to IRS spectra, in Mq. 

Note. — Errors listed are 90% confidence intervals from fitting to the post-shock density in column 2. Porosity, V, is percentage 
of grain occupied by vacuum. Three rows for compact grain case correspond to different values of shock compression ratio of 2.65, 4.8, 
and 6.8, as described in text. LMC abundances are assumed, so that Ue = I. Inn- Line strengths are for shocked ISM, and are reported as 
fractions of observed line strengths measured from RGS spectra by model with absorption described in text. For this remnant, observed 
Une strengths are: O Lya = 4.41 xlO"^^ gj-g^ cm"^ s~^; O Ka = 5.05 xlO-^^ gj-gg cm~^ s~^. Models predicting fractional line fluxes 
of greater than 1 are ruled out. M/3 is total amount of radiating dust, extrapolated from single slit position to entire remnant. M /j/Mq is 
ambient, pre-shock dust/gas ratio, where sputtering effects are accounted for. 



Table 3. Predicted ISM Oxygen Line Strengths for SNR BOS 19-69.0 



V 




no'' 






OLya 


OKa 




Mb'' 




Mc/Mg 









6.2i:i 

6 ?6.8 
"■^5.6 


2.34i:fl 

1.291:1? 


4.69|:^| xlO^° 
4.69|i| xlQio 
4.69ti| xlOi° 


Q Q9II.9 

5.47t:l§ 

3.86i:fl 


2.io?:# 

^■^"0.99 

0-82o:7o 


2.632:91 

^■^"1.30 

i-02j:^f 


15.6}!:} 

8.60?# 
6.06i:t| 


1.95?:|| xio- 
1.952:11 xio- 
1.95?:^i xio- 


Q 

— 

-3 
-3 


i.3|:f XlO" 

2-3?:^ xio- 

3.2i:^ xio- 


A 

-4 

-4 
-4 


25% 


Q 5IO.3 
='•-'8.6 


^ i:o3.89 
J--'°3.25 


7.186^:1^x1010 


23.2f^9f 


^ 0^4.29 

J-°J3.38 


4.i9i:?g 


23.9if:i 


93A^^jl XlO- 


-4 


3.9|:I xio- 


-5 


50% 


17I8.5 


fi 406.98 


1.29l:f^ XlO" 


74.6i:i 


7.80|:4| 


7.95f:il 


42.8i:6 


4.80t:i XlO" 


-4 


I.IJJt XlO- 


-5 



'^Hydrogen number density in cm~^. 
''Ionization timescale in cm~^ s. 
"^Emission measure (EM) in 10^^ cm~^. 
Swept gas mass in Mq. 
''Dust mass from fit to IRS spectra, in Mq. 

Note. — Same as Table 2, except that measured line strengths from RGS data are as follows: O Lya = 5.81 xlO~i^ ergs 
cm~2 s~i; O Ka = 3.91 x 10"^^ ergs cm~2 s~i. Models predicting fractional Une fluxes of greater than 1 are ruled out. 
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Fig. 1. — Top left: MIPS 24 yum image of SNR B0509-67.5, overlaid with regions of spectral 
extraction as described in text, where green region is "faint" region and yellow marks "bright" 
region. Top right: Star-subtracted H« image of 0509. Bottom left: MIPS 24 yum image of SNR 
B0519-69.0, with extraction slit shown in yellow. Bottom right: Star-subtracted Ha image of 0519. 
Spectral extraction regions in both remnants are approximately 10.5" x 45". FWHM of MIPS 24 
/im PSF is approximately 7". Ha images obtained with the 4-m Blanco telescope at NOAO/CTIO. 
In all images, north is up and east is to the left. 




15 20 25 30 

A (/zm) 



Fig. 2. — Top curve: Spitzer IRS spectrum of "bright" region of 0509. Bottom curve: spectrum of 
"faint" region. Note the differences in the slope of the spectra, in addition to the overall intensity. 
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Fig. 3. — 14-35 //m IRS spectrum of the "faint" region of 0509, overlaid with model fit. 
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Fig. 6.— Left: Ha image of B0509-67.5 from Figure [B convolved with the MIPS 24 /xm PSF. 
Right: MIPS 24 /im image, at native resolution. 



Fig. 7. — 24 yum image of 0509. Arrow shows both direction of slit placement and length of 
inferred brightness enhancement in that spatial direction. Diamonds are extraction centers for the 
parallel slit positions extracted in the text (note that arrow coincides with position 3). Intensity 
function plotted in Figure [8] corresponds to diamonds from left-to-right in this image. Green circle 
is MIPS 24/im PSF. 
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Arcseconds 



Fig. 8. — Intensity of background- subtracted integrated IRS flux as a function of distance from the 
left-most to the right- most diamond point in Figure |7l 




.5 IK 15 mi es: 



Wavelenglh (A) 

Fig. 9. — Top: XMM-Newton RGS spectrum of B0509-67.5, from 5-27 A, with model overlaid as 
described in text. Lines of interest are marked. High shock speeds and ejecta velocities lead to 
significant blending of Ne lines with numerous Fe L-shell lines. Data beyond 23 A was excluded 
from model fitting, and is shown only for ease of comparison with bottom panel. Bottom: Same, 
but for BOS 19-69.0. Slower speeds allow separation of Ne and Fe L-shell lines in this remnant. 
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0.0010 



0.0005 
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Fig. 10. — Best fit values for 0509 of post-shock density and dust mass, with errors, for various 
grain models. Errors are 90% confidence limits, as described in the text. 
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Fig. 1 1 . — Same as Figure[l0l but for 05 19 
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Fig. 12. — Evolution of the forward and reverse shocks with time for several values of the adiabatic 
index, 7, of the shocked medium. Top curves correspond to the forwa rd shock; bottom curves 
correspond to the reverse shock, t' and r' are normalized to values given in lPwarkadas & Chevaliei 
([l998). Solid line here corresponds to the solid line in Figure 2a of their paper. Note that decreasing 
7 has the effect of shrinking the post-shock volume. Effective 7 values of 5/3, 1.285, and 1.18 
correspond to effective r values of 2.65, 4.8, and 6.8, respectively. 



-38- 



SILICON DIOXIDE 




Energy (MeV) 



Fig. 13. — Stopping power of a proton in silicon dioxide (Si02) as a function of proton energy. 
Taken from the NIST PSTAR database (see text). 
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Fig. 14. — Ratio of proton heating rate of a si licate grain frorn calcul ations in this work to ap- 
proximations from lOraine & Salpeteil (| 19791) & lOwek & Wernen (| 198 1|) for an impinging particle 
energy of 100 keV, as a function of grain radius. 
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Fig. 15. — Comparison of spectra produced in dust heating model for protons of 100 keV. Solid 
line: spectrum assuming proton projected ranges calculated in Sectio n 7.5. Dashed line : spec- 
trum assuming analytical approximation to proton projected range from Praine & Salpeteii (|l979h . 
Electron heating (Te = 2 keV) is included in this model. Projected range is defined as the average 
length traveled by a particle into a material. 
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Fig. 16. — Ratio of spectra from Figure [T5l 
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A. Proton Energy Deposition Rates 

In young SNRs, and in 0509 in particular, proton temperatures are much higher than electron 
temperatures, and heating by protons cannot be neglected as it can in slower shocks where temper- 
atures have had time to equilibrate. Here we investigate whether simple analytic approximations to 
proton and alpha particle energy deposition rates are valid at ion energies approaching or exceeding 
100 keV. 



The stopping power (energy loss per unit path length) and projected range (average length 
traveled by a particle into a material) for both electrons and protons into a given grain are generally 
given by analytic expressions approximating experimental data. For electrons, thes e expr essions 
fit the experimental data to within 15% for the energy range of 20 eV to 1 MeV JPweklllQSTl) . 



For protons and alpha parti cles, the approxima te stopping power, based on experimen tal data from 



Andersen & Ziegleii (|1977|) and lZiegleii (|1977|) . is given by lDraine & Salpeted (I1979|) as 



?cp^ 



Rh = 3x IQ-^ p~\E/keV)g cm''^, Rue = O.GRh, (A1) 

where p is the mass density of the grain. The authors caution that this expression is valid only for 
< 100 keV. We raise here two questions: 1) Are these approximations appropriate for energies 
at or exceeding the 100 keV threshold (as may well be the case in very young SNRs with shock 
speeds exceeding 6000 km s~^)? 2) These approximations are based on data over 30 years old; are 
they still valid? 

To answer both of these questions, we turn to more recent laboratory data on the stopping 
power of protons and alpha particles in various materials, obtained from the online PSTAR and AS- 
TAR databases published and maintained by the National Institute of Standards and Technology^ 
In Figure [131 we show the stopping power of a proton in silicon dioxide as a function of energy. 
The stopping power turns over above ~ 100 keV. Since, for an extremely fast shock, there will 
be a significant population of protons with energies in this range, this turnover should be properly 
accounted for in calculating the projected range of the particle. 

Interstellar grains are more complicated than silicon dioxide (although it is possible that Si02 
is a minor component of ISM dust). Since the NIST databases only have information available for 
actual materials that can be measured in the lab, materials like "astronomical silicate (MgFeSi04)" 
are not available. The individual stopping powers for the constituent elements are available, how- 
ever, and can be combined to appro ximate the st opping power for a grain of arbitrary composition. 



To do this, we use "Bragg's Rule," (|Bragglll905|) . given by 



'Databases are available at |http : / /physics . nist . gov/PhysRef Pat a /Star /Text /contents . html| 
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S{A^B^)=m-S{A)+n-S{B), 



(A2) 



where S{AmBn) is the total stopping power of a molecule A^Bn, and S{A) and S{B) are the 
stopping powers of the individual con stituents. Since Mg is not included in the NIST database, we 
use the results of iFischer et alj (Il996|) . Thus, for MgFeSi04, one has 



S{MgFeSiOi) = S{Mg) + S{Fe) + S{Si) + 4 • S{0). 



(A3) 



The calculated projected range for protons based on this stopping power (similarly for graphite) 
are given by the following polynomial expressions in logarithmic- space, where E is the energy of 
the impinging particle in keV: 



Protons 



log {RH,siiicate) = 0.053 log^' (E) - 0.202 log' (E) + 1.21 log (E) - 5.68 (A4) 



log iRH,graphite) = 0.053 log=^ (E) - 0.158 log' (E) + 0.961 log (E) - 5.46 (A5) 
Alpha Particles 

log {RHe,smcate) = 3.78 X 10-^ log' (E) - 8.79 X 10^=^ log^ (E) - 4.25 x 10"' 
log^ (E) + 6.25 X 10"' log' (E) + 1.04 log (E) - 5.79 



log {RHe,graphite) = 3.33 X 10"^ log^ (E) - 6.13 X 10"^ log* (E) - 3.46 X 10- 

log^ (E) + 3.03 X 10"' log' (E) + 0.966 log (E) - 5.76. 



Ultimately, it is necessary to determine how much energy is deposited into the grain by an 
impinging particle, in order to calculate grain heating rates and determine whether these improved 
energy deposition rates matter. In Figure [H we show the ratio of the heating rate, H(a,T), 
calculated from the NIST energy de posit ion rates and that calcul ated from the analytical approxi- 
mations of iDraine & Salpeten (Il979() and lDwek & Wemeii (|198lh . for a proton temperature of 100 
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keV, which is relevant for the fast shocks, such as those seen in 0509. Figure [l5] shows a compar- 
ison of the spectra of grains heated behind such a shock, for both the modified energy deposition 
rates described above and the analytical approximations. Figure [16] shows the ratio of the spectra. 
The main difference is in the 18-40 fim range, but the maximum difference is only ~ 15%. The use 
of analytical approximations to proton heating rates should be valid for most cases. Nonetheless, 
we use the modified rates in reporting results in this work. 



A.l. Moving Grains 

Because heating by protons is non-negligible in the case of such fast shocks and high ion 
temperatures, it is necessary to account for relative gas-grain motions between dust grains and ions 
in the plasma. Grain heating rates for protons and alpha particles are increased by 35 and 112%, 
respectively, over the case where grains are at rest with respect to the plasma. Grain heating is 
still dominated by electrons, though, so the overall heating rate is increased by only ~ 25% for the 
moving grain case. When folded through the model, this results in a decrease of about 10% in the 
plasma densities necessary to fit IRS data. These numbers are reported throughout the text. We 
assume that the charge-to-mass ratio of the grains is sufficiently low that grains are not immediately 
affected by the passage of the shock. 



